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1. Introduction. 

Conformal field tlieory and integrable field theory in two dimensions are two subjects 
which attracted a lot of attention in the last years. The reason essentially lies in specific 
two-dimensional symmetries which lead to exact solutions of the quantum field dynamics. 
In conformal field theory solvability of a model is provided by a chiral algebra, which 
includes Virasoro algebra as subalgebra. Examples are W-algebras [00, parafermions 
1^ , Wess-Zumino model |^] , etc. When we go away from criticality all the symmetry is 
usually lost. However, as it was pointed out in 0, for certain classes of perturbations of a 
conformal field theory the resulting two-dimensional field theory possesses infinite number 
of mutually commutative conserved charges. The currents corresponding to these charges 
can be viewed as a deformation of some operators from (enveloping of) chiral algebra. The 
conserved charges result in a drastically simplified scattering theory - a general S matrix 
factorizes into product of two-particle S matrices [0 . Factorized scattering preserves the 
number of particles and the set of on-mass-shell momenta. This great simplification makes 
it possible to probe a vicinity of a conformal fixed point in detail [§]. Large variety of 
theories with factorized scattering was explicitly constructed |Tl[] |]12[ [|T3[ |T^ . Many 
of them are reductions fl^ at a special value of coupling constants of quantized classical 
field theories, e.g. Toda field theories. The factorization is typical for the scattering 
of solitons of the nonlinear classical field equations integrable by the inverse scattering 
method. The spectra of conserved charges in such a theories are essentially the same as 
in the reduced versions. Factorized scattering shows up in semiclassical limit as absence 
of backward scattering of lumps. Most of the considered models contain one or few stable 
particles and exhibit no resonances. 

In this paper we consider a model containing a bunch of unstable particles, besides a 
stable ones. From the point of view of conformal field theory the this model is a pertur- 
bation of SU{2)k Wess-Zumino- Witten model by a certain multiplet of primary operators. 
The perturbation depends in general on four arbitrary parameters. The layout of the 
paper follows. In section 2 we introduce the necessary notations and recall some basic 
facts about Wess-Zumino- Witten model. In section 3 we formulate the model, obtain re- 
currence formula for classical conserved currents and construct A^-kink solution for the 
model describing scattering, merge and decay of kinks and their bound states. Section 4 is 
devoted to the construction of quantum conserved charges. Besides the series of integrals 
of motion which survive semiclassical limit there is an additional set of conserved charges, 
not admitting semiclassical limit. A brief consideration of scattering amplitudes is given 
in section 5 
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2. Basics. 



In this section we recall basic facts about Wess-Zumino-Witten model and introduce 
necessary the notations. 

Consider two-dimensional Wess-Zumino-Witten (WZW) action fl^, |15[|,|1T9[] 



Sn.,^ = 1^/ ^^^^'^ {d^'gg-'d^gg-^) + kV [g] (2.1) 

where field g(t, x) takes value in any semisimple compact group G, and t, x - are coordinates 
of 2-dimensional Euclidean space {t is euclidean time) and k is some dimensionless coupling 
constant. Here we will deal mostly with G = SU{2). The ambiguous term F [g] is given 
by the functional 

r[^] = dtdxdre^^xTt [d^'-gr' d^-gg-' d^~gr') (2-2) 

where integration is over 3-dimensional half-space (t, x,r), r > 0. Boundary conditions 
g{t, x,oo) = 1, g{t,x,0) = g{t,x) define the functional F modulo 2ttN for some integer 
A^. With these boundary conditions functional integral is well-defined if /c is a positive 



integer number. It shown in 0|]2^the model (|2.1| ) is conformally invariant and therefore 



can be studied by methods of two-dimensional conformal field theory (CFT) [^. In fact 



this model possesses larger symmetry with respect to current algebra su{2)k x su{2)k |2C 



||22|| . Namely, the action (|2.1|) is invariant under transformations 

g ^ n{z)g{z,z)n{z) (2.3) 
where Q{z) and O(^) are arbitrary SU {2)-valned functions of light-cone variables 

z = {t + ix)/2, z = {t-ix)/2 (2.4) 
The field content, anomalous dimensions of the fields and equations for correlation 



function were studied in [22]. Partition functions and applications to string theory were 



discussed in |l2^] . The infinite dimensional symmetry ( p.3|) of the WZW model is generated 



by local currents 

J = -\dgg-\ J = -\g-^dg (2.5) 

which satisfy equations ||] 

dJ = dJ = (2.6) 
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The currents J{z) and J{z) obey the following operator product expansion (OPE) 



r{zi)j\z2) = ^^-^—^ + , . J\Z2) + regular terms (2.7) 

I \Z\ — Z2) [Zi — Z2) 

here tensors g"^ and /"^ are su{2) invariant metric and structure constants. They have 
components 

r = 1<1+- = 1, /o+- = 2, /r = -f-- = 1 (2.8) 

Let JF be space of mutually local fields of the theory. For any F{z, z) G we define an 
action of operators J^, (a = ±, and n = 0, ±1, ±2, ..) acting in this space 

J«F(^,^)= / ^,r[Q{C-zrF{zrz) (2.9) 
Jc-, 27rz 

As it follows from ( p.7|) operators satisfy commutation relations of su{2)k Kac- Moody 
algebra 

[Jn^J^l = fc^Jn+m ^Q^^^n+mfi (2-10) 

Since the theory has conformal invariance energy-momentum tensor is traceless and has 
only two independent components 

(2.11) 

Field T{z) satisfy OPE 

/ X / X c 1 2T(z2) dT(z2) 
T{zi)T{z2) = --, TT + 7 + 7 ^ + regular terms 2.12 

2 [Zi - Z2)^ (^1 - Z2y [Zi - Z2) 

Again as in the case with Kac-Moody algebra we can define action of generators L^, n = 
0, ±1, ±2, .. in the space JF 

L„F(z,z)= / ^TiOiC -zr+'F{z,z) (2.13) 
Jc, 27rz 

Operators Ln form Virasoro algebra 

[L^.Lm] = {n-m)Ln+m + —in^-n)5n+jn,o (2.14) 
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with central charge 

c = ?,k/{k + 2) 
and commute with Kac-Moody generators as 



[J, 



(2.15) 



(2.16) 



In terms of Virasoro generators Ln can be expressed as 

+00 



qab 



{k + 2) ^ 



b 

n—m 



(2.17) 



m= — 00 



where symbol : (■ ■ ■) : means normal ordering, i.e. operators J" with n < to be put to the 
left. Operators J" and L„ related to the fields J°'{z), T{z) can be defined in a similar way . 
There are invariant fields [|4[, [|2l$J;^-^^, jjj = 0, 1, 2, .., k/2, m,m = -j, -j + l, ..,i-l,i 
in the space which obey equations 



J«$H'S = 0, n > 
J^^^l^'l = [{j -m){j + m + 1)]'/^ 



'0 m,m 
"^0 m,m 



+ l,m 



(2.18) 



(7 7") 

and the "bar" -counterpart of the above equations. Invariant fields ^jii in have dimensions 
(A,.A5), 

A, = ^ (2.19) 

If i = j we will denote the corresponding field as The space can be represented 

as a sum 



(2.20) 



where 



is a space spanned by fields of the form 



jai 7-02 Tfls T-bi 7-&2 jbi 



^2.21) 



In semiclassical limit k ^ oo fields become an SU{2) matrix g^^\z, z) in represen- 

tation of spin j while fields ( |2.21| ) become differential polinomial in matrix elements of 
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3. Perturbation of WZW model. Semiclassical analysis. 

Now we turn to perturbation of WZW model. As a perturbation we talce operator O 

O = ^ j dtdx {Tr{gAg-^B) - Tr (AS)) (3.1) 

In this section we make rotation t —>■ —it to Minkowski space. The perturbed action 
functional become 

A = S^,^ + O (3.2) 

here A and B are diagonal matrices 

A = Diag(ai,a2,a3) , S = Diag (6i, 62, 63) (3.3) 

and we take g in real spin j = 1 representation. Equivalently we can say that g takes values 
in SO{3) instead of SU{2). If we introduce coordinates 1/;, on SO{2>) group manifold 

(cos ip sin ip \ / cos tp sin Lp\ / 1 ^ \ 

-sinip cosip 1 cos^? sin^? (3.4) 

1 j \ — sin (f cos cp J \0 — sin 1} cos 1!} J 

the action functional A takes the form 

^= dx~^dx~ (^8+4) d-'ip + d+if d-if + d+^) d-'d — 2 sin ipd-{) dj^^! — U{'tp,ip^^}^ 

(3.5) 

where the potential U is of extremely complicated form 

U = sin^ i/j + a2/32sin^'i? + aa/^a sin^ 99 + 

00 00 00 

a2/9i sin '0 sin 1} — asPi sin 1/; sin ip — a2l3s sin sin + (3.6) 
a2/3i sin?/' sin(/? sint? (2 cos '0 COS"!? + sin sin sin i?) , 

x^ = {t ± x)/2, = 4//c and coupling constants defined as 

CKl = ai - 02, ^2 = ^2 - ^3, "3 = fll - 0^3, 

(3.7) 

Pi =bi- 62, /32 = &2 - &3, P3 = h - h 



The dynamical model described by the above action was first considered in [53 



For general and f3i the perturbation completely breaks su{2)k x su{2)k symmetry 
of the model. Another feature of the model is completely broken P-parity , i.e symmetry 
with respect to refiection of space coordinate x —>■ —x. WZW can be made invariant under 
change of sign of x if we simultaneously perform transformation g g~^ ■ In the perturbed 
model we do not have this option because of the specific form of the operator O. 



3.1. Lax pair and classical integrals of motion (IM). 



In [^a Lax pair representation for (|3.5| ) was found and it was shown that the model 
can be solved using inverse scattering method p^. For our purposes we will use another 
Lax pair which is more simple and suitable for analysis. 

Equations of motion for the perturbed model (|3.2|) 
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are compatibility condition for an auxiliary liner problem 



d+J = [gAg-\B] 

k ^ (3.8) 

J = --rd-gg 



£_^' =(d- - d-gg-^ - iX'^ kB^^ 
C+-^ = (^d+ - iXKgAg-^^<i' = 



(3.9) 



where A is a spectral parameter. It is remarkable enough that both the equation ( |3.8| ) 
and the operators C± (but without the spectral parameter A) appeared also in a context 



of = 2 supersymmetric models The Lax pair (|3.9| ) is a guaranty for existence of 



infinite number of conserved charges 



-Q(^) = 



(3.10) 



in the model (|3.2|) (at least on classical level). In Q^"^ s refers to the spin of the charge and 



a distinguish between charges of the same spin. Using the Lax pair we find the recursion 

(s) 

formula allowing to construct the densities Qa 

= (jK^^A (3.11) 

V / aa 

where matrix K^^^ defined recursively 

Ki'J ^ 

V / mn 

Sl+S2=S 

with initial conditions 

= 0, = (3.13) 
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A few first densities are 

Q['^ =Tt{JJ), 
Qi''> =IV(SJJ), 

^(3) _ _ \ ^ JgpJpqJqsJsa JgpJpa Jgp Jpa 1 41 

r. JaqJqpd—Jpa d—Japd—Jpa 

^ ih,-ha){h^-haY ^ (bp -bar 
p,q^a ^ ' ^ f pj^a ^ 

(3) 

Only two of Qa are linearly independent. Further analisys shows that allowed set of 
spins is s = ±1, ±3, ±5, ±7, ±9, ... and for any value of s there are two linearly inde- 
pendent conserved charges. Densities with negative s can be obtained from the positive 
ones replacing J — > jg~^d-^.g and bm CLm- In the next section we construct quantum 
conserved charges corresponding to ( |3.14| ) . 

3.2. Kinks. 

Now we are going to discuss soliton content of the model. First of all we specify 
boundary conditions. The natural ones B^rc q ^ QzLoo x ^ zbcxo, with QzLoo 

being some 

constant matrices. They can be determined from the condition of vanishing of the right 
hand side of ( |3.8[ ) as x ^ ±oo. There are four allowed vacua 

go = Diagonal (1, 1, 1) 
gi = Diagonal (—1, —1, 1) 

(3.15) 

(72 = Diagonal (1, —1, —1) 
gs = Diagonal (-1, 1, -1) 

All these vacua have same energy. Therefore we expect that there are kinks interpolating 
between them. Let for a moment f = and ^? = in ■ Then we get familiar 

sine-Gordon (sG) action for the field ?/'• It is well known that sG theory has kinks ||28|| , 



and the corresponding solutions interpolate in our model between go and gi. The same 
happens if we set another couple of fields to zero. However, there is significant difference 
between sG kinks and kinks in the model we are considering. Because tti (5*0(3)) = Z2 
kinks must have Z2 charge. For example lets take kink g{t, x) interpolating between go 
and gi. At given moment t, g{x) is a path on S'0(3) manifold connecting identity element 



with Qi. But gi is nothing but rotation by angle ±7r and kinks with different Z2 charge 
correspond to different choice of the sign. On the contrary, sG kinks have Z charge. 

Probably the easiest way to construct kink solutions is to use Riemann problem |2^] for 
the Lax pair ( |3.9|) . In this approach solitons are related to \E''s in (|3.9|) having simple 
analytical properties in A. Namely, 



'^{x ,a;''~|A) = Il{x , a;''~|A) ■ exp (zA ^kBx + iXnAx^) 



(3.16) 



where 11 is a meromorphic function of A. If we normalize 11 as n(x~, a;^|0) = 1, then an 
example of such a solution for 11 relevant for scattering of kinks and having poles in 
A-plane can be written as 

detM("^)(A) 



[n(A)]„, 



detM 



(3.17) 



where matrices M^"'')(A) and M are (A^ + 1) x (A^ + 1) and N x N correspondingly 



M 



7i + 7j 



(3.18) 



and 



(b) 



A el") 



M("^)(A) 



2— \ 



'1 



A+i7„ 



M 



(3.19) 



here ■ji are positive numbers (they are related to rapidities 9 of kinks) and (e^, ej) is a 
scalar product of vectors 



e-a = exp(7„KAx"'~ — 'j^^kBx ) ■ 



(3.20) 



here Cn are some real 3-dimensional vectors independent of z and z. 
The corresponding solution for g is 

g{x~,X^) = 5fa ■ n(x", X+|oo) ■ fiffc 



(3.21) 



and the current J is determined by expansion of 11 near A = 

[n(x-,x+|A)]_ = 5_ + ^A ^^l^jf^^]- + O(A^) 



(3.22) 



The simplest solutions are kinks Kq^]{6) with topological charge e interpolating be- 
tween go and Qa- They correspond to the case = 1 and one of components of the vector 
c IS zero in ( 13:21] ) . For example, for i^^^'V) 

sinh^i (-l)'^l 



cosh 01 cosh < 

_sinh^ I (3-23) 

cosh 01 cosh 01 


here 

<Pi ^ («i7 + /3i7"') 2: + I (ai7 - /3i7"') t + Si (3.24) 

ei = 0, 1 is Z2 topological charge of the kink and 61 is an arbitrary real constant. We 
see that only field if^ is exited. There are similar solutions for kinks Kq'^2\9) and Kq'^^\9) 
when only fields ^? or are excited. Using the energy-momentum tensor we find that 
kinks i^ilV) and ones obtained from them by multiplication by g^s ( see (|3.15|) ) have 
two-dimensional momenta 

Pj"^^ = (M„cosh^,,M„sinh^,) (3.25) 

where masses of kinks 

Mr, = -{kaM^/^ (3.26) 

and rapidities 9n are related to 7 

= -ln7 + ^ln(/3^/a^) (3.27) 

now we can rewrite phase ( |3.24| ) as 

27rMi 

In the following we will assume that 



(xcosh^i - tsinh^i) + Si (3.28) 



ai > 02 > 03, 61 > 62 > 63 (3.29) 

In general we have 3A^ free parameters in (|3.21| ) . Restricting components of 
Ci by condition c\^^ ■ c\^^ = 0, z = 1,2, ..,A^ we obtain 2N parametric solution de- 
scribing scattering of kinks. We give an example of scattering of two different kinks 
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4f{92)Kil'\9^) {N = 1, cf) 



0, c 



(1) 



0). From (gj) 



and ( |3.21| ) we find, for example 



tan-i/; 
simp 
tant? 



-(-1)- 

(-1) 
-(-1)- 



2(A - l)e^^+'f>- 



(3.30) 



1 -j- 62-^1 — A^e^''^^ 



32((^l+(^2) 



here 



A = tanh 



27rM„ 



A; 



6*1 — 92 



(x cosh 6*7^ 
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a2/3i, 
t sinh^^) + Sn 



(3.31) 
(3.32) 



In full form the corresponding solution is given in Appendix A. 

What are the rest of parameters responsible for? In general solution (|3.21|) de- 
scribes scattering, merge and decay of solitons. The simplest example is decay process 



(£l+e2) 
03 



-^02^'' (^2) -^23^^ (^1) (-^ — 1' components of Ci are nonzero) 
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with phases 



_2( — e<t>l+24>2 
l_|_e2<^.2 _|_g2(</)j^+c^2) l-|_g2<^2 -|-e2(<Al+'#>2) 
_2(-l)"=le*l+^*2 l_e2,A2_|_e2(<#.i+</.2) 
T+^2^?2+^20?7T?27 i^e2'''2 -|-e2Wl+'*2) 



-2(-l)^ 



-2(-l)"2e'^ 



2(-l) 



<!l+'!2p<#>l+<#>2 



2(-l) 



<!2 p<P2 



^e2<#>2 _|_g2((#.j^+<^2) l-|-e2'''2 -|-e-^('^l+'^2) 



l^e2'''2 -|-e2(*l+'^2) 
]^_g2<;)2 _g2{<;>]^+(^2) 
l-j_e2<#>2 -|-e2(<#'i+<#'2) 



(3.33) 



(3.34) 



Indeed, suppose that 0:37 — /337 ^ = 0, 027 — /327 ^ > 0, and consider worldline 
(pi 4'2 = const. Then if t ^ —00 

/ i_e2(,^i+,^2) -2(-ir 
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lJ^e'^(4'i + 't>2) 





2(-l) 



<!l+<^2 p'l>l+4>2 



l_|_e2(</'i+<^'2) 





1 





"T+^2T7T+^2r 



l_e2{<;>i+<;)2) 



l_|_e2{,#,i+02) 

i.e. at t = —00 we have resting kink Kq'^^~^'^^\ Now consider limit t 
(pi = const. In this limit we get 

9 — ' 



-2(-l)^ie'^i 



v 







-2(-l)'^le'^l 







\ 




(3.35) 

+00. We first keep 
(3.36) 



1/ 
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that is i^23^'' kink propagating in the positive x direction. Taking hmit with (f)2 = const 
we obtain -^"02^^ kink propagating in the negative x direction. 



/I 













2(-l)"2e'»2 



. ^ 

-2(-l)"2e'^2 
l+e2</'2 

/ 



(3.37) 



That instabihty of the heaviest kink Kq^ can be also predicted from the mass spectrum of 
sohtons ( ^ ) 

Ms > Ml + M2 (3.38) 



The solution 



can also be obtained 



where equality achieved when A = B in (| 
from ( 13.301 ) in a limit 71 ^72- 

Now we have a very peculiar situation: the model has infinite number of local con- 
served charges (|3.10|) and at the same time there are unstable particles. If all particles were 



stable the standard argument [prelate the presence of infinite number of local conserved 
charges to factorization of a scattering matrix into two-particle S matrices. In quantum 
theory unstable particles are usually excluded from "in" and "out" states, and the only 
trace of them are resonance poles in S'-matrix of the theory. From this point of view it is 
crucial to find quantum counterpart of charges (|3.10|) to discuss quantum integrability of 
( p72|) , S'-matrix, etc. We address this problem in the next section. 

As in sG model there are bound states of kinks. Solutions describing bound states, 
scattering processes of bound states and kinks can be obtained from (|3.21|) by replacement 
of some poles 7^ by pair of poles 7^ ± iSn- 

3.3. so{2) reduction of the model. 

One can note that when one of q:„ or vanish, the model regains right or left so{2) 
symmetry. If, for example, as and Ps vanish simultaneously so{2) x so{2) symmetry is 
regained in the model. It means that there is a propagating massless field. It is possible 



to factor out this mode as in ||3T[], [^and get another model described by the functional 

d^ud,,u 



S = / d^zi ^ + 2A ((uu) 

/ \ 1 — uu 



uu 



(3.39) 



This model is closely related to Z„ parafermionic models perturbed by second thermal op- 



erator £2 [p3 |. 
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which after change of variables 



w = exp (^y^ sin , u = exp (^-y^ sin , (3.40) 

takes the form 

S = d^zl^d^edf.e + tan^ (9/2) d^ad^^a + A(cos26' -1)) (3.41) 
This model is integrable and possesses conserved charges of spin s = ±1, ±3, ±5, .... To 



describe the densities of these charges we introduce new fields u and 7 which satisfy 
equations 

duj = da — TT, Buj = da 



2cos2|' 2cos2f 

1 n (3.42) 

1 K cost/ 



97 = da —o", 97 = da- 



2cos2f' ' 2cos2f 

The above equations are compatible if equations of motion of ( |3.41|) are being imposed. 
Using ui and 7 we construct currents 

j± = (d9 ± idatsiu^j e^'^ 
j± = [de ± i^atan- e^'^ 



\ 2 
For A = equations of motion for (|3.41| ) are equivalent to 

dj^ = dj^ = (3.44) 

Local integrals of motion (IM) for ( |3.41| ) turn out to be differential polynomials in and 
j^. The first ones, s = ±1, ±3 are 

Q^'^ = rr 

9 3.45 

Q^'^ = (i+r) + 2dj+dj- 



2 
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4. Quantum Integrals of Motion. 

In this section we will construct IMs proving quantum integrability of the model 
discussed in section 3. 

The central object of quantum theory are correlation functions of some operators On 

< Olizi,Zl)---OniZn,Zn) >A= j V(j) Ol{zi, Zl) ■ ■ ■ On{Zn, Zn) BXp {-A[(p]) (4.1) 

where fields (p ^-re some "basic" fields, operators are some local functionals of ^'s 
and their derivatives, and A [(p] is an action functional. In this section we are back to 
euclidean two dimensional space. In our case A is sum of WZW action ( pJ] ) S^zw, 
which is conformally invariant and represent fixed point of renormalization group, and 
perturbation ( ^.1| ) . We will assume that the structure of space of fields in the perturbed 
model is the same as in WZW model (|2?20D . To identify operator O as an operator in 
WZW model we use (13. 61) to obtain 



(4.2) 



= (9(2) ^ 0(0) 

where vectors ry"^, fjrh and constant rj are given by 

?7 = -TiA ■ TiB - Tr [AB) 
3 



(4.3) 



(2) 

For > 6 operators ^m\m ^'^^ relevant and we will show soon that the action is well defined 
for /c > 10 (see discussion of " resonance" condition below) . After expansion in O^"^^ to 
the first order we get 

< Cl(zi, Zl)--- On{Zn, Zn) >A = < C'l(^l, Zl) ■ ■ ■ On{Zn, Zn) >wzw 6^ - 

< Ol{zi, Zl)--- On{Zn, Zn) >wzw + ■ ■ ■ 

(4.4) 

Here < ■ ■ ■ >wzw means that average is taken with respect to WZW action and V is 
a volume of the base space. Note that the perturbation theory around conformal fixed 
point is IR divergent and to make it finite we fix V finite but large. Here we concentrate 
on equations of motion which are relations between local operators and their derivatives. 
Prom this point of view IR divergences are irrelevant. 
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Suppose we want to calculate d derivative of some holomorphic (at fixed point) oper- 
ator, say Q{z), which is local with respect to $^^1(2;, z). After taking derivative of (ET 



and neglecting usual contact terms (coming from the first term in (|4.4| ) ) we get ||3 

dQ{z,z) = ^rnVrn I dC (C, ^) Q{z) (4.5) 

If the residue of OPE 

Res[r]rnfirh^^^,fh{z,z)Q{w)] = dR{w,w) (4.6) 

V / z='w 

for some operator R, then we have constructed quantum conserved current 

dQ{z,z) = dR{w,w) (4.7) 

In general d derivative of a local operator Q in perturbed model schematically take a form 

BQ = {r]ff)Ri + {mfR2 + ■ (4.8) 

where Rn are some local fields. Because the structure of the space of fields in the perturbed 
WZW model was assumed to be the same as in unperturbed model, we conclude that the 
series in the right hand side of (|4.8|) terminates at some term. Moreover, all fields R^ must 
have anomalous dimensions coinciding with the WZW ones. If the anomalous dimension 
of some field Rn matches the one in WZW it must be added to the right hand side of 
( [4.5|) . The "resonance" condition can be derived from ( f4.8|) by matching dimensions of 
operators on both sides of the equation 

1 - n (1 - As) = A, (4.9) 

for some n > 1 and j. The solutions to the above equation corresponds to additional term 
of order [rjff)^ in ( [4.8| ) . The only "resonance" cases are 

/c = 6 j = 0, n = 4 and j = 1, n = 3 
/c = 7, J = 0, n = 3 

(4.10) 

/c = 8, j = 1, n = 2 
A; = 10, j = 0, n = 2 

For that values of k one need to consider next to the leading order of the conformal 
perturbation expansion ( [4.4| ) to derive equations of motion for local operators. Therefore 
we will restrict our attention to the case /c > 10 when the action is well defined and will 
work within the first order of the conformal perturbation theory. 
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4-1- Coserved currents in U (^su{2)f^^. 

Here we will find conserved currents which are composites of Kac-Moody currents 

Following 0]we define define new operators Dn^m-, n = 0, ±1, ±2, .., m = —2, 2 by 
its action in current module spanned by fields Jn\Jnl---Jnl ' ^ 

D^,nQ {z, -z) = i dC (C, z) {C- zT Q{z) (4.11) 

where notation $m, [z, z) for the operator rjrh^rn^m {z^ z) was introduced. Using ( |2.1(j| ) 
one can easily prove the commutations relations 

[j^p,Dra,n] = mDm,p+n (4-12) 
[j~,Dm,n] = [{j + m){j + Dm-l,p+n 

These operators satisfy simple relations 

VmDmfiQiz) = -dQ{z,z), 

1 (4-13) 

Now it is straightforward to find d derivatives of any field from the current module - 
apply rjmDmfi and then move all D^s all the way to the right where they transforms into 
derivatives of ^rn\m.- Sugavara construction ( |2.17|) allow to express the result entirely in 
terms of Kac-Moody generators J^. Solving after that (|4.7|) we find quantum integrals of 
motion for spin s = 1,3. For example 

; (4.14) 

The expressions for first two nontrivial quantum integrals of motion of spin s = 3 are given 
in Appendix B. 

What we found is in perfect agreement with the classical result (|3.14|) . As before 
there are 2 conserved charges of spin s = 1 and 2 charges of spin s = 3. Charges with 



m • 



s = —1 and s = —3 can be obtained from ( |4.14| ) and replacing —J^ and rjm fj. 

It is reasonable to make 

Conjecture. For general k the model ( p.2|) possesses infinite number of mutually 
commutative local integrals of motion of spins s = 0, ±1,±2,... and for any s there are 
two independent IM. 

Note that the derivation does not feel if k is integer or not. 
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4-2. Additional set of conserved charges for integer k. 



In this section for even k we derive additional set of integrals of motion. 

The spectrum of anomalous dimensions of WZW theory (|2.19| ) contains Ak = k/A. 
It is possible to show (see Appendix C) that the symmetry algebra of WZW model can be 
extended by new currents 



= and ^^{-z) ^ 



(4.15) 



For k = 4n, for example, these fields contribute to the partition function of SO{3)4n WZW 
model, [EBI. We are going to show that certain linear combinations of give rise to 

new set of IM. 

To construct new IM's we need to consider OPE, (14.61) , 



C(2, -, - - 2) ry^ry^ ^— -^^ + ^——^ 



Z2, Z2\ 



+ 



(4.16) 



Here C(2, 2, 2 ~ 2) is a structure constant of OPE (we will not need its exact value), 
(3l^^ri and are some coefficients that can be determined comparing transformation 

properties of the left and right hand sides of ( [4.16| ) under action of su(2)k. To find „ 



and Vi introduce following p4|] new variables x,x and generating functions 
J{x,z) = J-{z) + 2xJ^{z) - x^J+{z) 

1/2 ^, , (4.17) 



^'^^'^\x, x\z, z) 



J.J 



^rn+j -m+j ^Ud) 



z,z) 



m=-j,m=-j 

where C^+^ = {2j)\/{j + m)!(j — m)\ are binomial coefficients. The singular part of the 
necessary OPE's (|2.7|) take a form 

J, .J, X AX1-X2Y ^Xi-X2 J, . {Xi- X2Y dJ{x2,Z2) 

J{Xi,Zi) J{X2,Z2) = -k- — -2 J{X2,Z2) - 



{zi - Z2y 

J{xi,Zi) $^-^'^Hx2,X2|^2,^2) = 



Zl - Z2 



Zl - Z2 



dx-: 



2j- — + ^- — ^^^'^'ix2,X2\z2,Z2) 



Zl - Z2 Zl - Z2 dx2 



(4.18) 



As in (|2.9|) we define action of modes J"(x) 

1 



J'^{x)F (x, x\z, z) 



(4.19) 
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where e+ = Eq = —S- = 1. OPE of primary fields can be written in a form 



yc I' J2 J 



3,J 



(4.20) 



(X2, X2\Z2, Z2 



where 



$(^'^^)] {xi,Xi\zi,Zi) = n]^^^^izi2,Xi2)^^j^{zi2,Xi2)'^^'~'\x2,X2\z2,Z2) (4.21) 



and Zi2 = z\ — Z2, x\2 = xi — X2, etc. Operators T^-^j^ defined by series 

oo 

7^3^,,,(^l2,a;l2) = J] r\{x2)r\{x2) ■ ■ ■ Jl\{xr,) Rn{ai, ..,ar,) 

n=0 ai=0,±l 

where Rn{ai, .., an) are some differential operators 

Rn{ai,..,an) = ^rp(ai, .., an|xi2) J 



(4.22) 



(4.23) 



Operators Rn are analogs of coefficients and in ( |4.16| ) . There is, of course, 

similar expression for TZi -. . In the following to simplify notations we will often omit all 
"bar" arguments and indices as we will work in the holomorphic sector. 

Operators i?n(ai, ••, cin) can be obtained, as we already mentioned, recursively from the 
requirement that both sides in (|4.20|) transform in same way under su{2)k. In particular 
to get the recursion we apply J"(x) with n = and n = 1 to both sides of (|4.20| ) . Action 
with Jq{x) relate different r^'s within same R^ and action with J"(x) allow to relate rp's 
for two successive R^. The operator Rq was found in [^. For (ji, 725^) = (2, |, I — 2) it 
is identity operator 

Ro = 1 (4.24) 

It is possible to find the next operator i?i(a). The only nonzero coefficients rp{a) for the 
case (ii, J2,i) = (2, |, | - 2) are 

2{k-6) 



ro(0) = 
roi-) = 
riH = 



k-4 ' 
4 



k-4' 
2 

k-A 
17 



(4.25) 



Now, finally, we introduce functions 

2 



fix) = fp^'^" 

p=-2 
k/2 

g{x) = J2 9px''/^+'' (4.26) 

p=-k/2 
k/2-2 

h{x) = hpx'^^^-^+P 

p=-k/2+2 



to define operators 



$(2)^^) _ j dfi^^\x,x)f{x)^^^\x,z), 

^'(z) = j dfx'^^\x,x)g{x)^^^\x,z) , (4.27) 

$(|-2)(;2) = f c^^(t-2)(x,x)/l(x)$(t-2)(a;,^) 



where (i// is a measure 

d^i^^\x,x) = - (4.28) 

We have arrived at the problem: when the residue Res [$^2)(^-|^^v[f(22)] ^ is a total 
derivative 



Res 



The equation ( |4.29| ) is equivalent to the system of linear equations 

dn^^Hx,x)g{x) (^{k-6){a + l)p{x)x'' +p'{x)x''+^ - 2p{x)x''+^ ^|-2+m 



2 A; + 2 



j dii^^-^\x,x)h{x) (2:"+'^ - ^(A;-4)(a + l)x"^ x^" 



for a ={0,±1}, m = {-| + 2,---,|-2} 

(4.30) 

where 

p{x) = /2 - hx + /ox2 - /_ix3 + /_2x4 (4.31) 

This system has 3/c — 9 equations for 2k — 2 variables Qp, hp, i.e. for k > 8 the system is 
overdetermined. We have found solutions of this system for f±i = and /c = 6, 7, 8, ■ ■ • , 24. 



18 



Such a choice of /p's corresponds to the perturbation {'^^ . From analysis of obtained 
solutions one can draw 

Conjecture, (i) For any k = An, 4?! + 1, 4n + 3 and k >11 the model ( |3.2[ ) possesses 
two additional IM which are certain linear combinations of currents '^m- (ii) For any 
/c = 4n + 2 > 14 the model (|3.2|) possesses three additional IM which are certain linear 
combinations of currents ^m- 



Despite that the system (|4.30| ) has also solutions for "resonance" cases /c = 6, 7, 8, 10, 
we can not say anything about conservation of the corresponding charges in these cases. 
For first nonresonance cases k = 11, 12, 13, 14 we give solutions for Qp^s in Appendix D. 



5. Decay of particles. 

In section 3 we have constructed the solution describing decay process for the kink with 
mass M3. The mass of that kink is always larger or equal (if A = B) then sum of masses of 
other, lighter, two kinks ( J3.38D . It is natural to conjecture that while in the corresponding 
quantum theory the heavy kinks are being excluded from "in" and "out" states they result 
in resonance poles in S'-matrix. To justify this point we calculated semiclassical limit of 
^-matrix for the scattering i^ii'^(^i) k[1^\02) ^ K'^t\^2) i^23'^(^i), (PD 

siZ^l ^ out < K^2'\e'2)K^,^\e\)\Ki\^\e,)Kti\e2) 

( 1 S{e'i - ei)d{9'2 - 62) exp (i (^KoiKi3 - ^K,n - ^i^is) j 

where Ag is action functional A evaluated on the corresponding kink solution ( |3.23|) , 
( p.30|) . For any solution g of ( |3.8|) integrals on the left hand side in (|5.1| ) can be done 
with the use of formula for a variation of A 

^Ma] = -^TrJ dxdtd^ {i9-'Sg)ig-'d^g)) (5.2) 



One can obtain 



, mk k r /l + e^i2-eo-iC 
exp / aCln 



27r Jn \ e^i2-^o -|- 



(5.3) 



here = i In ((ai/32)/(a2/3i)), ^12 = 61-62. From (U) follows that ^J^J^^^g^ (6*) has a 
line of essential singularities in the physical strip < Im6 < tt at Re6 = Oq. We expect 
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that this singular hne corresponds to series of resonance poles which become dense and 
are "glued" together as /c — oo. 

We can also study scattering of "basic" particles tjj, "d, (f which enter the action func- 
tional A. From (|3.5|) we read their masses 

= 2(fc-iai/3i)i/2 

= 2(A;-ia2/32)'/' (5.4) 

As in the kink case we obtain 

m^f, > + ms (5-5) 
This is an indication of instability of the ip particle. Indeed, to the leading order we get 



out 



< ^{p^mP2)W{pi) >in=i(^] k-\ail32 - a2l3i)5^^\pi - p2 - Ps) (5.6) 



In the rest frame of particle ijj particle propagate to the right and "d propagate to the left. 
Due to broken P-parity in the theory the amplitude for decay process with interchanged 
t/j and ^? is zero. It is possible to check that to the leading order amplitudes for process 
■0 ^ is zero. This illustrates integr ability in the tree approximation - nontrivial IMs 
prohibit such kind of processes. The amplitude for scattering t/j and "i? is 

SM<^12) = out < i^{e\)^{9'2mo2mei) = 

here sign ± distinguish between -i/; particle propagating to the left or to the right in the 
center-of-mass frame. Due to existence of higher IMs in the theory reflection part of the 
amplitude is absent in (|5.7|) . The shift of rapidity in (|5]^) agrees with the one in 
semiclassical amplitude for kinks (|5.3|) . 
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6. Conclusions. 



In the paper we studied quantum integrability of SU{2)k WZW model perturbed by 
multiplet of primary fields with isospin j = 2. The model contains four free parameters. 
The spectrum of the theory contains kinks and their bound states. Kinks interpolate 
between different vacua in the theory and has Z2 charge. Some of kinks are unstable 
and we obtained solution describing decay of the heaviest kink into two lighter kinks. 
Quantum integrability was proven by explicit construction of quantum IMs. It appears 
that there are two sets of IMs. The first one - IMs belonging to universal enveloping algebra 
U (^su{2)i^ . Integrals of motion from this set admit semiclassical limit k ^ 00 and there 
is recurrent formula for IMs in this limit. The other set are integrals of motion built from 
chiral components of primary field of isospin j = k/2 (conformal dimension Afe = j). 
Obviously these IMs do not admit semiclassical limit. We also briefly discussed scattering 
amplitudes for kinks and fundamental particles. We argued that scattering amplitudes for 
stable kinks and particles possess resonance poles and calculated to the leading order in 
1/k decay amplitude for the heaviest particle. The work on the construction of the exact 
kink S-matix of the theory is in progress. 
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8. Appendix A. 

Full expression for the two kink solution is given by 



911 

913 

922 

931 
933 



1 _ e^<Pi 


+ A^e^'^^ 




1 + 62-^1 

2(-l) 


_|. A2e2<^2 4- e2(<^i+02) ' 


1 + e2'^i 
1 - e^-^i 


+ A2e2'^2 
- A2e2<^2 


_^ e2(<^i+<^2) ' 

_^ g2(<^l+<^2) 


1 + e2'^i + A2e2<^2 
-2(-l)^i+^2(l - 


_|. g2(<^l+</>2) ' 

A)e'^2+<^i 


1 + e2<^i 
1 + e^'t'i 


+ A2e2<^2 
- A^e^'^2 


_|_ g2(<;ii+<^2) ' 

_ g2(<^i+<^2) 


1 + e2<^i 


+ A2e2<^2 


_|_ g2(<^i+<^2) 



912 
921 
923 
932 



2(_i)6ie'^i(l - Ae2<^2) 

1 _|_ g20i _|. A2e2'^2 + e2(<^i+<^2) 

-2(-l)^ie'^i(l + Ae^-^^) 
1 _|_ e2cj,i _^ A2e2'?^>2 + e2(<^i+'A2) 

-2(-l)^2(g2^1 _ ^)g02 

I _j_ ^24)1 _j_ ^2g2(j!>2 _|_ g2(<^i+<^2) 

-2(-l)^2(^A + e2<^i)e'^2 

1 _|_ g2<;!.i _|_ A2e2<^2 -I- e2{ct>i+<t>2) 



(8.1) 
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(9.1) 



9. Appendix B. 

Here we give expressions for first nontrivial IM's of spin s = 3 in universal enveloping 
algebra U ^sw(2)^^ 

Qf^ =(QV-2J-iJ-iJ-iJ-i + l2V6r]-2VoJIiJIiJ-iJ-i + 
20r]^2V2JliJliJ-iJ-i + IQV-2V2JZ1J-1J-1J-1 + 
4(977^ + 2r]-2V2)J-iJ-iJ-iJ-i + 12V6voV2J-iJ-iJ^iJ^i + 
6r]ptiJtiJtiJti + 12\/6 (8 - k) r]-2VoJl2J~iJ-i + 
16(fc - 10) ?7_2?72 JI2J-1J-1 + 8(/c + 2) r]-2V2 J^JIiJ^i + 
6^6 (12 - k) r]or]2 J^Jl^J^i + 6\/6 (/c - 10) 7]-2Vo JI2J-2 + 
{9{e - 12k + 16) ril + 2{e -8k -4) r]_2m) J-2J-2 + 
8(A;2 - la/c + 24) r/_2r/2 JZ2J-2) " 1 



Q2''' = - ^V-2 <^-i<^-i<^-i<^-i + '^2V6ri-2VoJ-iJ-iJ-iJ-i - 

4(677^ + 577_2?72) JIiJIiJ+iJli + 8(3r/2 - 277.2772) JZ^J\j\Jli + 

4(377^ - 277-2772) J\J\J\J\ + 12V6riomJ-iJ-iJ-iJ-i - 

677! J^iJ^iJ+iJ^i - 72V6r]_2VoJl2J-iJ-i + 

(12(/c - 16)77^ - 2(/c + 2)?7_2772)J°2-^ri-/^i + 

8(/c - 10) (377^ - 277.2772) JZ2J-1J-1 + 36v^77o772 J^J-i-^-i + 

(6(6 + 10k - k'^) ril - S{k'^ - 13k + 24) 77.2772) Jl2<^-2 + 

(S{k^ -8k- 40) 7/^ - 2{k^ -8k -A) 77-2772) J°2<^-2 + 

3V6{k - 10) (fc - 2) 77_2r/o ^13-^-2 + 

3\/6(A; - 10)(A; - 2) 770772 J^-^i) • 1 



(9.2) 
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10. Appendix C. 

In this appendix we prove that for integer k current algebra sw(2)^ can be extended 
by operators ( [4.15| ) . 

In section 4.2 we introduced generating functions J{x,z) and ^^^'^\x,x\z,z). Using 
Knizhnik-Zamolodchikov equation and null-field JZi{x) z) one can derive [^equa- 

tion for four-point conformal block, 

G =< ^^^\xuZi)^^^\x2,Z2)^^^\x3,Z3)^^^\x4,Z4) > (10.1) 

The solution of the equation takes a simple form 

Gll{x,z) = \x32X14ix - z)] (232^14^(^-1)) (10-2) 

where x = {x 12X34) /{x 14X32) , z = {Z12Z34)/ {Z14Z32). To check associativity of the al- 
gebra generated by ^^^\x,z) we fuse $'s in different ways and compare results. For 
example, we can fuse pairs $^^^(xi,2i), ^^i\x2j Z2) and (^3, ^3), ^^i\x4^ Z4), or 
^^i\xi, zi), ^^i\x3^ Z3) and ^^i\x2, Z2), ^^i\x4, Z4). To make functions Gl1{x,z) and 
6^24(1 — x,l — z) coincide one should impose "equal-time" commutation relations 

^^i\xi, zi)^'^^\x2, Z2) ^ e^^'^i\x2, Z2)^^^\xi, zi) , (10.3) 

With this prescription algebra generated by operators become associative. 
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11. Appendix D. 

Here solutions for functions g{x), ( |4.26| ) are given; cko,±i, ck and (3 are free parameters. 
= 11 

92^ = ^/o(80/oV-2/2-7/o^- 320/!2/|)/3 
f| = 3/2 « 

91 = ^/-2(16/oV-2/2-/o' + 64/!2/|)/3 

9^ = -5/0/2^ a 

^3 = -^/o/!2(/o+24/-2/2)/3 

^1 = /|(5/o + 8/,/2)« (11.1) 
= /^2(5/o + 8/,/2)/3 

^-1 = -^/o/l(/o+24/-2/2)« 



^-1 = -5/0 /-2/^ 

= Y^/2(16/oV-2/2-/o' + 64/!2/|)« 

9-^ = ^/o(80/2/-2/2-7/o^- 320/!2/|)« 
= 12 

^76 = Y^fl (8 /o' /3 + 27 /_2 /2 /9 + 15 /2 a) 
^4 = ^/-2/|/o(8/3 + 15a) 
92 = f-2 fl 

^7o = /-2/l/o« (11-2) 

^7-2 = /^2 fl /3 

^7-4 = ^/^2/2/0 (8/3 + 15 a) 

^7-6 = ^/-2 (8 /o /3 + 27 /_2 /2 /3 + 15 /o' a) 
fi'-5 = 9-z = 9-1 = 9i = 93 = 95 = 

24 



k = 13 

9^ = ^(3/o'-40/oV-2/2 + 192/oW|-512/VI)/? 

^1 = Y^/o/-2(/o'-16/oV-2/2 + 192/VI)/3 
= 4/o/| « 

^1 = ^/^2(/o -48/oV-2/2 + 64/!2/|)/? 

^1 = -/|(7/| + 8/_2/2)« 
^1 = 2/o/i2(/2+8/_2/2)/3 
= 2/o/|(/o' + 8/_2/2)« 

- -/l2(7/o' + 8/-2/2)/? 

= ^/|(/o - 48/0V-2/2 + 64/VI) « 

= 4/0/^2/3 

= Y^/o/2(/o' - I6/0V-2/2 + 192/!2/|) a 

^-^ = -^/^2/3 

= -^(3/o^-40/oV-2/2 + 192/oW|-512/VI)« 
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k=U 

97 = Y^/l (5 /o «i + 86 /-2/0/2 «i + 35 /0V2 «-i + 77 /_2/| a-i) 
96= Y^h "0 

^5 = T^/-2/2 (/o «i + 15 /-2/2 «i + 7/0/2 a-i) 
loo 

5 /• /•2 
^4 = - — /O/2 «0 

^3 = 4/-2/l(/o«i + 7/2a-i) 
15 

£^2 = /2 (/o + /-2/2) ao 

91 = /-2/I «i 

3 (11-4) 

90 = -^ fo ifo + 6 /-2/2) ao 

^-i = /-2/l«-i 

^-2 = /-2 (/o + /-2/2) ao 

9-3 = ^/-2/l (/o a-1 + 7/-2 ai) 

5 r2 ^ 

S'-4 = -YY/-2/oao 

^-5 = Y^/-2/2 (/o «-l + 15 /-2/2 CX-1 + 7 /-2/0 «l) 

5 ^3 
9-6 = J-2 «0 

9-7 = Y5^/-2 (5 /o «-i + 86 /-2/0 /2 «-i + 35 /_2/o' ai + 77 /i^^ ai) 
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